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ABSTRACT 

O ■ 

The mid- infrared local luminosity function is evolved with redshift to fit the 
O I spectrum of the cosmic infrared background (CIRB) at A > 5 /im and the galaxy 

^ I counts from various surveys at mid-infrared, far-infrared and submillimeter wave- 

lengths. A variety of evolutionary models provide satisfactory fits to the CIRB 
and the number counts. The degeneracy in the range of models cannot be bro- 
^ I ken by current observations. However, the different evolutionary models yield 

\ approximately the same comoving number density of infrared luminous galax- 

ies as a function of redshift. Since the spectrum of the cosmic background at 
A > 200 iim is quite sensitive to the evolution at high redshift, i.e. z > 1, all 
models that fit the counts require a flattening at 2; ~ 0.8 to avoid overproducing 
the CIRB. About 80% of the 140 fim CIRB is produced at < z < 1.5 while 
only about 30% of the 850 fim background is produced within the same redshift 
range. The nature of the evolution is then translated into a measure of the dust 
enshrouded star formation rate density as a function of redshift and compared 
with estimates from rest-frame optical/ultraviolet surveys. The dust obscured 
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star formation rate density appears to peak at z = 0.8 ± 0.1, much sooner than 
previously thought, with a value of 0. 25^q\^ MqI/t'^ Mpc^^ , in excellent agree- 
ment with a value derived from extinction correction to NICMOS observations, 
and remains almost constant up to 2; ~ 2. At least 70% of this star formation 
takes place in infrared luminous galaxies with > 10^^ Lq. The long wave- 
length observations that constrain our evolutionary models do not strongly trace 
the evolution at z > 2 and a drop off in the dust enshrouded star formation rate 
density is consistent with both the CIRB spectrum and the number counts. How- 
ever, a comparison with the infrared luminosity function derived from extinction 
corrected rest- frame optical/ultraviolet observations of the Lyman- break galaxy 
(LBG) population at 2; ~ 3 suggests that the almost flat comoving star formation 
rate density seen between redshifts 0.8 and 2, extends up to a redshift of ^ 4. 

Subject headings: cosmology: diffuse radiation — infrared: galaxies — galaxies: 
evolution 



1. Introduction 

The Extragalactic Background Light (EBL) in the infrared, also referred to as the Cos- 
mic Infra-Red Background (CIRB), is a record of the emission, absorption and re-radiation 
of photons integrated over cosmic history. It provides a valuable constraint on theories of 
galaxy formation and evolution. The EBL at near-infrared wavelengths is due to redshifted 
radiation from stars. At mid-infrared (MIR) wavelengths, the background is due to redshifted 
emission from dust that consists of the polycyclic aromatic hydrocarbon (PAH) features and 
very small grains transiently heated to T ~ 200 K in individual galaxies. At far-infrared 
(FIR) wavelengths, the dominant contributor is thought to be cold dust (T ~ 20 K) that 
is heated by the ambient interstellar radiation field in galaxies. The recent detection of 
this background at 2.2 /xm, 3.5 /im, 140 fira and 240 /xm using DIRBE/COBE data, and 
in the 125—2000 fim range using FIRAS/COBE measurements by various groups (Puget et 
al. 1996; Dwek & Arendt 1998; Fixscn ct al. 1998; Hauscr ct al. 1998; Schlegel et al. 1998; 
Lagache et al. 1999; Gorjian, Wright & Chary 2000; Wright k Reese 2000; Wright 2001), 
has indicated that the intensity of the optical/near- infrared background is roughly equal to 
that of the far-infrared background. This implies that about 50% of the integrated rest-frame 
optical/ultraviolet (UV) emission from stars and other objects is thermally reprocessed by 
dust and radiated at mid- and far-infrared wavelengths. Thus, star formation rates that are 
derived from rest-frame optical/UV luminosities of galaxies are a lower limit to the true star 
formation rate (e.g. Madau, Pozzetti, & Dickinson 1998; Meurer et al. 1999; Steidel et al. 
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1999; Yan et al. 1999). 

The first good evidence of this came from the IRAS sky survey which revealed a new- 
population of galaxies with Lir — L(8 — 1000 jim) > IO^^Lq (see review by Sanders & 
Mirabel 1996). Those with Ljr > lO^^L© were classified as ultraluminous infrared galaxies 
(ULIGs) while galaxies with 10^^ > Lm > lO^^L© were classified as luminous infrared 
galaxies (LIGs) ^. These objects exhibited the largest known star formation rates of all 
local galaxies but had ~90% of the bolomctric luminosity being emitted in the far-infrared 
(40—500 fim) indicating that dust reprocessing is a significant parameter that needs to be 
considered in estimates of star formation in certain galaxies (e.g. Soifer et al. 1986). However, 
in the local universe, the integrated bolomctric luminosity density of "normal'' optically 
selected galaxies is Lboi = 4 x 10® Lq Mpcr^ while that of infrared luminous galaxies is 
~ 8 X 10^ Lq Mpc'^ i.e. 50 times less (Soifer et al. 1987). This seemed to indicate that the 
contribution from LIGs and ULIGs is sufficiently small that they need to be considered only 
as extreme cases. 

Spectroscopic follow-up of the faint IRAS population which covered a relatively small 
redshift range {z < 0.27) indicated that infrared luminous galaxies were more numerous 
in the past than they are today and may have a significantly larger contribution to the 
integrated luminosity density than inferred from observations of the local universe (Kim & 
Sanders 1998). Deeper observations that trace the far- infrared luminosity of galaxies to 
high redshift are difficult since cirrus and confusion noise rapidly begins to dominate. 

The ISOCAM guaranteed time extragalactic surveys in conjunction with the ELAIS 

survey and observations of the lensing cluster Abell 2390 covered a range of flux densities 
between 50 yuJy and 50 mJy at 15 /im (Altieri et al. 1999; Elbaz et al. 1999; Serjeant et 
al. 2000). The differential counts resulting from these surveys revealed that the counts of 
galaxies increase quite rapidly as 5*^7'^ at brighter ffux levels {S^ > 0.4 mJy) and then flatten 
out as S~^'^ at fainter levels. The observed mid-infrared counts are an order of magnitude 
higher than expected if the local mid-infrared luminosity function were not evolving with 
redshift. This rapid increase in mid-infrared luminous galaxies has been modeled as a (1 -|- 
zY^^ luminosity evolution in the 15 /im local luminosity function by Xu (2000), and as 
a combination of number density and luminosity evolution by Franccschini et al. (2001). 
This evolution is much stronger than observed in the UV by Cowie, Songaila & Barger 
(1999) who flnd that the comoving UV luminosity density evolves as (1 + z^''^ instead of 
the (1 -|- 2;)^-^='=°-^^ initially proposed by Lilly et al. (1996). Furthermore, observations of 



•^Previously, the term LIG was used for all objects with L/h > IO^^Lq. We use "infrared luminous 
galaxies" when referring to both LIGs and ULIGs collectively. 
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galaxies in the local universe have shown that the mid-infrared and infrared luminosities 
are well correlated (Section 2). The mid-infrared luminosity of ~ 70% of the sources seen 
in the ISOCAM surveys, particularly in the Hubble Deep Field-North and flanking fields 
(HDFN-I-FF) translates to an infrared luminosity greater than 10^^ Lq implying that the 
majority of them are LIGs and ULIGs (Elbaz et al. 2001). At 2; ~ 0.8, the mid-infrared 
luminosity density derived from the ISOCAM 15 /xm sources is ~7xl0^ Lq Mpc~^ while 
the 8 — 1000 /im luminosity density adopting the mid to far-infrared correlation seen in the 
local universe is ~5xl0^ Lq Mpc~^. In comparison, at z ~0, LIGs and ULIGs contribute 
10^ Lq Mpc~^ to the 12 jim. and 15 jjtm luminosity density and 7.8x10^ Lq Mpc~^ to the 
infrared luminosity density, as derived from the local luminosity function (LLF) of Soifer 
et al. (1987); Fang et al. (1998); Xu et al. (1998). This indicates an increase by a factor 
of about 60 between z ~0 and z ~0.8, providing further evidence for an evolution in the 
infrared luminosity function (IRLF) with redshift. 

Similar deep surveys have been conducted at 850 fim using the SCUBA instrument 
on the JCMT (Hughes et al. 1998; Barger, Cowie, & Sanders 1999; Blain et al. 1999a; 
Bales et al. 2000). The large beam size (14" FWHM) and the negative /c-correction in this 
wavelength regime makes identification of the optical counterparts and thereby the redshift 
distribution of the sources very difficult. High resolution radio interferometric observations 
and the use of 450 //m/850 /xm flux ratios have helped somewhat in locahzing the sources 
and constraining the redshifts (Hughes et al. 1998; Barger, Cowie, & Richards 2000). These 
have placed the bright (5",^ > 6 mJy) submillimeter sources at 2; ~ 1 — 3, which has been 
confirmed by the more extensive survey of Chapman et al. (2001). The imphcation of this is 
that most, if not all the submillimeter sources are extreme ULIGs with star formation rates 
of 10^ — 10'^ Mq yr^^. Furthermore, the star formation rate density due to ULIGs must have 
increased by about 2 orders of magnitude between z ~0 and z ^ 1 — 3. 

Many of the LIGs and ULIGs in the local universe show morphological signatures of 
interaction, and >50% of the optical counterparts of ISOCAM HDFN galaxies show evidence 
for interactions (Mann et al. 1997). Surveys at visible wavelengths show a redshift evolution 
of the merger fraction, defined as the fraction of close pairs of galaxies, as ~ (1 -|- z)^ (e.g. 
Le Fevre et al. 2000). Thus, if mergers were indeed a tracer of LIGs and ULIGs, this would 
again suggest that the bright end of the IRLF is evolving strongly. However, the faint end of 
the IRLF is very poorly constrained at 2; ~ 1 since none of the long wavelength surveys are 
sensitive enough to detect galaxies with Lm < 10^^ Lq eit z > 0.5. Meurer et al. (1999) has 
shown that the FIR to UV flux ratio is closely related to the UV slope for normal starbursts 
but that the relationship breaks down for ULIGs (Meurer et al. 2000). This indicates that 
the visible/near-infrared counts can potentially place constraints on the evolution of the faint 
end of the infrared luminosity function but we postpone this discussion to the future. 
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In this paper, we combine data from a variety of published surveys of nearby galaxies to 
determine the correlation, if any, between the luminosities at various mid- and far-infrared 
wavelengths. Wc use these correlations to generate smoothly varying spectral energy dis- 
tributions for galaxies as a function of luminosity class. We assess the need for luminosity 
and density evolution in the 15 /im luminosity function of Xu et al. (1998) and therefore 
the 60 iim luminosity function of Soifer et al. (1987) based on fits to the ISOCAM 15 /im, 
ISOPHOT 90 /xm, ISOPHOT 170 //m and SCUBA 850 fim galaxy counts, as well as the 
spectrum of the CIRB at A > 5 /im. The evolution of the mid-infrared luminosity function 
is then translated to an estimate of the dust enshrouded star formation rate density (SFR) 
as a function of redshift and compared with SFR values derived from optical/near-infrared 
surveys. We adopt a Ho=75 km s~^ Mpc~^, Qm — 0.3, = 0.7 cosmology throughout this 
paper unless otherwise explicitly stated. 

2. Luminosity Correlations in the Infrared and Template Spectral Energy 

Distributions 

It can be shown that the 12 /xm and far-infrared luminosities of galaxies in the IRAS 
Bright Galaxy Sample (BGS) cannot be accurately derived from their 5— band luminosities 
(Soifer et al. 1987)"^. The peak to peak scatter in the Lb/Liji ratio for a fixed Ljfi is about a 
factor of 20. However, as mentioned earlier, the FIR to UV flux ratio has been shown to be 
closely related to the UV slope for normal starbursts (Meurer et al. 1999). This relationship 
breaks down for the ULIGs (Meurer et al. 2000). The phenomenon can be qualitatively 
explained by the fact that the UV emission arises from stars which are relatively unobscured 
to the observer. Regions of star formation with large optical depth i.e. H II regions, could 
exist in the central regions of galaxies where almost all the UV hght is reprocessed to the 
MIR and FIR. Thus, the regions of FIR and UV emission would be unrelated, especially for 
ULIGs. The best observational evidence for this explanation can be seen in the Antennae 
galaxy (Mirabel et al. 1998) where about half the 15 fim emission seen by ISOCAM arises 
from regions that are inconspicuous at visible wavelengths. The breakdown in the FIR to 
UV slope correlation for ULIGs is problematic for determination of the true star formation 
rate from optical/UV surveys, since submiUimeter observations using the SCUBA instrument 
indicate that ULIGs might have a larger contribution to the star formation rate density at 
high redshift. This suggests that it will be difficult to determine the true SFR by applying 
an accurate extinction correction to the optical/UV determined value. 



"^The Zwicky magnitudes, m^, in the BGS were converted to B-luminosities using ms = rUz — 0.14 and 
a -B— band zero point of 4260 Jy. 
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Fig. 1. — Plots showing the relative accuracy of tracing infrared luminosities (8—1000 /xm) of 
IRAS BGS (Soifer et al. 1987) objects from the 5— band (0.44 /xm) and 12 /xm luminosities. 
The lower plots show the ratio between true infrared luminosity {Ljji) and predicted infrared 
luminosity (< Lm >) derived from the 5— band or 12 /xm luminosity using a first order 
polynomial fit to data points with Ljr > 10^° Lq. If the infrared luminosity of all galaxies 
could be predicted precisely from their S— band or 12 ^m. luminosity, then all points in the 
lower plot would lie in a horizontal line with Lir/ < Lm >—l. The lower plots also show 
the la uncertainty in the prediction of the infrared luminosities. 
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Fig. 2. — Plots showing the relative accuracy of tracing infrared luminosities (8—1000 /xm) 
from the 6.7 /im and 15 luminosities. The 6.7 /xm luminosities are from ISOCAM guar- 
anteed time surveys (filled circles; Chanial et al. 2001) and ISOPHOT 7 fim spectroscopy of 
ULIGs (triangles) by Rigopoulou et al. (1999). Note that the Rigopoulou et al. (1999) values 
have been modified as described in the text. Asterisks are the 6.7 //m luminosities for the 
starburst dominated ULIG sample of Tran et al. (2001). Some of the extreme ULIGs might 
have a significant AGN contribution which could result in a deviation from the 6.7 //m-to-IR 
correlation derived for starbursts by decreasing the IR luminosity for a given 6.7 luminosity, 
as for the two brightest galaxies in our sample. The 15 /xm luminosities are from the Chanial 
et al. (2001) (filled circles) and the Aussel et al. (2000) (pluses) sample of galaxies. The 
lower plots are similar to those in Figure 1. 
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Since the short wavelength starhght and dust emission are not closely related, an esti- 
mate of the dust enshrouded star formation rate can only be derived from other tracers such 
as the mid-infrared and far-infrared luminosities or using the radio to far-infrared correlation 
shown by Condon (1992). In the mid-infrared regime, the spectra of galaxies exhibit broad 
emission features at 6.2, 7.7, 8.6, 11.3 and 12.7 fim which are probably from polycyclic aro- 
matic hydrocarbons (see review by Puget & Leger 1989). These features and their associated 
continuum dominate the emission at mid-infrared wavelengths shortward of 10 /im. There is, 
in addition, a continuum from very small grains (VSGs) of size <10 nm which dominates the 
emission above ~10 iim (Desert et al. 1990; Laurent et al. 2000) except for quiescent star 
forming galaxies. The VSGs get transiently heated to temperatures of ~200 K by the ambi- 
ent optical/UV continuum which is proportional to the star formation activity. In addition, 
mid-infrared measurements do not need large extinction corrections since the extinction at 
mid-infrared wavelengths is only about 1% that at visible wavelengths (Mathis 1990). The 
radio wavelengths on the other hand are dominated by free-free emission from HII regions 
and synchrotron emission from supernova remnants. Although radio observations are almost 
confusion hmited at an 8.5 GHz sensitivity of 9 //Jy obtained over the HDF (Richards et 
al. 1998), we find that they are typically as sensitive as the ISOCAM 15 //m observations in 
that they can typically probe galaxies with Lm ~ 10^^"^ at 2; ~ 1. Since the ISOCAM 
15 //m observations provide the primary constraint on cvohition models at ^ < 1.2, we adopt 
as a starting point the local 15 /xm luminosity function described in Xu et al. (1998) and Xu 
(2000). For the rest of this paper, we will use the convention defined in Sanders & Mirabel 
(1996): 

LiR = L(8 - 1000 fxm) = 1.8 x 10"^^ x 10^^ [13.48Li2 + 5.I6L25 + 2.58L60 + Lioo] (1) 

LpjR = L(40 - 500 fim) = 1.6 x 1.26 x 10"^^ x 10^^ [2.58^60 + i^ioo] (2) 

In the above, the symbol Lx is defined as L^{X /im) in units of L© Hz~^. Lm and LpiR are 
in Lq. 

To use the mid-infrared luminosity function as a tracer of the dust enshrouded star 
formation rate, we first need to define a calibration scale. Figures 1 and 2 illustrate the 
accuracy with which the infrared luminosity of galaxies can be derived from their mid-infrared 
luminosities. Figure 1 is based on ~300 galaxies from the IRAS BGS while Figure 2 is based 
on published ISOCAM and ISOPHOT observations of IRAS galaxies. The data points with 
LiR > 10^° Lq are fit by a first order polynomial shown in Equations 4-6. This is shown as 
a solid black line in the upper panels of the figure. Objects with Ljr < 10^° Lq are not used 
for the fits. At these low luminosities, the fraction of the bolometric luminosity emitted in 
the infrared is <50% i.e. visible starhght that is not obscured by dust is dominating the 
radiated energy of the galaxy. On the other hand a significant fraction of the more luminous 
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objects show disturbed morphologies, suggesting interactions with other galaxies that would 
result in gas-rich systems with star formation in highly obscured regions. The lower plots in 
both figures show the scatter in the ratio of the infrared luminosity as derived from IRAS 
data for these galaxies to the infrared luminosity derived from the polynomial fits. Also 
shown is the la uncertainty in the derived infrared luminosity calculated as the range within 
which 68% of the galaxies lie. The lowest luminosity objects {Lm < 10^ Lq) have been 
rejected in the lower plots. 

The 6.7 fim luminosities were derived from ISO observations of ~90 nearby starburst 
dominated galaxies. 44 spiral and starburst galaxies had photometry from ISOCAM (Cha- 
nial et al. 2001; Roussel et al. 2001; Laurent et al. 2000), 8 ULIGs had spectra from 
ISOCAM circular variable filters (CVF) observations (Tran et al. 2001) while 37 ULIGs had 
ISOPHOT mid-infrared spectra (Rigopoulou et al. 1999). Rigopoulou et al. (1999) obtained 
mid-infrared spectra with ISOPHOT of about 60 ULIGs and about 15 low luminosity star- 
bursts and normal galaxies to study the emission features from the PAHs. Of the 60 ULIGs, 
about 45 had the 7.7 /xm PAH feature detected with good signal to noise. However, the 
calibration and performance of the instrument is not very well determined. To assess the 
quality of the data set, we compared the ISOPHOT observations on 5 ULIGs and 4 low 
luminosity starbursts/normal galaxies to the Chanial et al. (2001) ISOCAM LW2 observa- 
tions of the same galaxies^. We find that for the ULIGs, the ratio of 7.7 /xm line-|-continuum 
fiux density as pubhshed in Rigopoulou et al. to the ISOCAM 6.7 //m fiux density lies in 
the range 1.5—3.0. In comparison, for the starbursts/normal galaxies, the ratio of 7.7 /im 
hne-|-continuum to the ISOCAM 6.7 /xm fiux density falls in the range 0.2—1.6, a factor of 8. 
So, for assessing the correlation between the mid and far-infrared luminosities, we consider 
the ISOCAM data on bright IRAS galaxies as well as the ULIG sample of Rigopoulou et al. 
(1999), dividing the line+continuum fiux value published in the latter by 2.4 and assigning a 
peak to peak error bar of a factor of 2. This is consistent with the range of 1.5—2.7 that we 
find for the 7.7 /xm lineH-continuum to the 6.7 /xm fiux density ratio in the ISOCAM CVF 
observations of Tran et al. (2001). 

The 15 fim and infrared luminosities of 120 IRAS galaxies were taken from the sample of 
Chanial et al. (2001) and the survey performed in the North Ecliptic Pole Region (NEPR) 
by Aussel et al. (2000). The NEPR sample of galaxies only has IRAS 60 iim luminosities 
available and we have converted these to a far-infrared luminosity based on a 60 /im to 
far-infrared correlation derived by combining the IRAS BGS and the IRAS PSCz catalog 
of Saunders et al. (2000). The far-infrared luminosity is typically about 83% of the total 



^The LW2 filter is broad enough to include the 7.7 /xm PAH feature but is centered at 6.75 /xm. 
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infrared luminosity and we have applied this conversion to be consistent with the other plots. 

Clearly, the 6.7 /xm, 12 /xm and 15 /xm luminosities of galaxies trace the infrared luminos- 
ity much better than the S— band luminosity. The 15 /xm-to-IR and 12 //m-to-IR correlations 
were determined from a first order polynomial fit to the data with Lir > 10^° L©. Since the 
6.7 //m-to-IR correlation is based on different data sets from CVF and broadband photometry 
with ISOCAM and spectroscopy using ISOPHOT, the polynomial fit for the 6.7 ^m-to-IR 
correlation was determined by applying a A;— correction from 6.7 fim to 15 fim based on ISO- 
CAM observations of nearby galaxies (Figure 3) and then using the 15 //m-to-IR correlation. 
The mid-infrared to IR correlations show a similar scatter around the correlation line which 
is about a factor of 5 better than the optical-to-IR correlation. 

These data sets tentatively illustrate the potential of using the mid-infrared as a tracer 
of dust enshrouded star formation and a more homogeneous and comprehensive survey of 
nearby galaxies, as will be undertaken by SIRTF, will be required to either strengthen or 
reject this correlation. 

Kennicutt (1998) has transformed the infrared luminosity of young < 10*^ yr starburst 
galaxies to a star formation rate. If we adopt the correlations shown in the previous fig- 
ures, we can translate the mid-infrared luminosity of galaxies with Lm > 10^° Lq to an 
approximate estimate of the dust enshrouded SFR (p') using the formulae: 



where all values are in Lq. The la values have been estimated by calculating the range of 
values within which 68% of galaxies have their observed infrared luminosities. 

As mentioned earlier, the main observational constraints on models that trace the red- 
shift evolution of the IRLF are: 

1. differential counts from various surveys at mid-infrared, far- infrared and submillimeter 
wavelengths and 

2. the spectrum of the CIRB at A > 5 /xm. 




(3) 
(4) 
(5) 
(6) 




To use these constraints, it is necessary to know the luminosity at different wavelengths 
for galaxies in each luminosity bin of the infrared luminosity function. This motivates the 
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generation of template spectra for objects of different luminosity classes. It is useful to 
note that many evolutionary models which have already been developed either use a mid- 
infrared template which is a very poor representation of the true PAH emission features (e.g. 
Rowan-Robinson 2000; Pearson 2001; Sadat, Guidcrdoni & Silk 2001) or neglect the PAH 
features altogether (Malkan & Stecker 1998, 2001). It is trivial to show that this makes 
a critical difference in the quality of the fits to the ISOCAM mid-infrared number counts 
and therefore affects the evolution parameters particularly at z < 1. So, it is important to 
have template SEDs that reproduce the observed trend in the luminosity of local galaxies at 
different wavelengths. 

Using MIR, FIR and submilhmeter data from ISOCAM, IRAS and SCUBA observa- 
tions of nearby galaxies, we fit the observed trend between different mid- and far-infrared 
luminosities as shown in Figure 3. The top two panels in the figure show the ISOCAM 
observations at 6.7 /im and 15 /im for ~50 IRAS galaxies which are described in Chanial 
et al. (2001). The solid lines in the panels show two first order polynomial fits, one for 
galaxies with 15 fim luminosity less than 2 x 10^ Lq and another for more luminous galax- 
ies. This is because the ratio between the mid-infrared luminosities chang function 
of 15 nm luminosity, possibly due to enhanced emission from the VSG component (Laurent 
et al. 2000). The luminosity break corresponds to Li^ ~ 2 x 10^° Lq which is similar to 
the luminosity cutoff used for deriving Equations 4-6. The panel showing the 15 /xm to 
60 //m trend consists of data described in Figure 2. The 60 //m to FIR correlation is for the 
IRAS BGS and PSCz galaxies while the panel showing the IR to FIR correlation is only for 
the IRAS BGS galaxies. The last panel shows SCUBA submilhmeter data on ~100 IRAS 
galaxies (Dunne et al. 2000). In addition, the correlation between IRAS 25 fim and 100 //m 
luminosities for galaxies in the BGS was also determined. The solid lines for the four lower 
panels utilize only a single first order polynomial fit to all the data points. Also shown in the 
panels as red triangles arc the luminosities at the corresponding wavelength for the different 
templates that were generated as described below. 

Template spectral energy distributions (SEDs) were generated between 0.1 and 1000 fim 
to reproduce the observed trend between mid-infrared and far-infrared luminosities. To 
generate these templates, we used the basic Silva et al. (1998) models to reproduce the 
ultraviolet-submiUimeter SED of 4 prototypical galaxies: Arp220, NGC6090, M82 and M51. 
These correspond to objects of four different luminosity classes - ULIGs, LIGs, "starbursts" 
(SB) and "normal galaxies" respectively. ISOCAM CVF observations between 3 and 18 //m 
of these galaxies provided new data on the relative strength of the mid-infrared features 
and continuum (Charmandaris et al. 1999; Laurent et al. 2000; Forster-Schreiber et al. 
2001; Roussel et al. 2001). The mid-infrared region of the modelled spectra were then 
replaced with the ISOCAM observations. In addition, corrections were made for the 17.9 /im 
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Fig. 3. — Plot showing the data (asterisks) at different wavelengths from IRAS, ISOCAM 
and SCUBA surveys (Soifer et al. 1987; Aussel et al. 2000; Dunne et al. 2000; Saunders 
et al. 2000; Chanial et al. 2001). The lines are the best fit polynomial of order 1. The red 
triangles are the corresponding values from our template spectral energy distributions which 
were generated as described in the text. 
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silicate feature based on observations by Smith et al. (1989). The four template spectra were 
checked to ensure that the IRAS observed values of these four galaxies were reproduced. 
We then partitioned the four templates into a mid-infrared (4—20 /im) and far-infrared 
(20—1000 ^m) component and interpolated between the four to generate a range of mid- 
infrared and far-infrared sample templates of intermediate luminosity. An additional set 
of far-infrared templates provided by Dale et al. (2000) were added to the ensemble of 
far-infrared templates to span a wider range of spectral shapes. 

For each luminosity bin of the 15 fim luminosity function, the luminosity at the following 
wavelengths: 6.7, 12, 25, 60, 100, 850 fim was predicted based on the polynomial fits to 
the data shown in Figure 3. Of the ~ 100 mid-infrared sample templates generated as 
described above, the mid-infrared template that best fits the predicted 6.7, 12 and 15 /im 
luminosities was selected. Similarly, the far-infrared template that best fits the predicted 
25, 60, 100 and 850 //m luminosities was selected. The luminosity of the templates at 
the corresponding wavelengths was determined by integrating over the filter curves of the 
instruments. Our goal was only to generate SEDs that reproduce the observed trend in 
luminosities at different wavelengths. Selecting a variety of sample templates provided better 
fits to the predicted luminosities than by just interpolating between the four SEDs generated 
by the Silva et al. (1998) models. The best fitting mid-infrared and far-infrared templates 
were then merged together to provide the final template SED for each luminosity bin. The 
red triangles in Figure 3 are the luminosities at the corresponding wavelengths from the final 
merged template SEDs. The S— band luminosity of galaxies in the IRAS BGS shown in 
Figure 1 was also used to constrain the optical/near- infrared spectral energy distribution of 
galaxies but as stated before, we have not constrained the UV slope of the template SEDs. 
The absence of a good correlation between the fi— band and IR luminosities implies that 
the optical/near-infrared part of our SEDs is highly uncertain. This is not a major problem 
since we arc only analyzing the dust emission in this paper. The templates for three objects 
with infrared luminosities of 10^°, 10^^ and 10^^ Lq along with the predicted luminosities at 
different wavelengths based on the correlations in Figure 3 are shown in Figure 4. 



3. Evolution of the 15 //m and Far-Infrared Local Luminosity Function 

Since our intention is to use the different mid- and far-infrared observational constraints 
to estimate the evolution of the dust enshrouded star formation rate with redshift, we use 
the 15 /xm LLF as a tracer of dust emission in the local universe. Xu et al. (1998) and Xu 
(2000) derived a 15 iim LLF based on a correlation between ISOCAM mid-infrared and IRAS 
mid- and far-infrared data. In addition, estimates of the 12 /im LLF have been made by 
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Fig. 4. — Template spectral energy distribution for objects of three different infrared lu- 
minosities along with the predicted luminosities at different wavelengths (diamonds). The 
luminosities correspond to L/^ = 10^^, 10^^ and 10^° Lq illustrating the SED of ULIGs, LIGs 
and starbursts respectively. The lower plot shows the same templates, normalized at 0.44 \im. 
(S— band) to show the evolution of the spectrum as a function of infrared luminosity. No 
correction for the UV slope has been made. 
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Rush, Malkan & Spinoglio (1993) and Fang et al. (1998). The different 12 //m and 15 //m 
luminosity functions arc shown in Figure 5. Also shown is the predicted 60 //m LLF derived 
from the mid-infrared LLF using a mid-infrared to 60 /xm conversion from the polynomial fit 
to the observations of Aussel et al. (2000) and Chanial et al. (2001) described earlier. All 
these are in good agreement with each other since they were essentially derived from IRAS 
observations of nearby galaxies. 

Evolution of the luminosity function with respect to redshift can be expressed as: 



where ^(L, z) is the number density of galaxies as a function of luminosity L and redshift 
z. The n(z) term represents evolution in the number density of galaxies while the 
term represents luminosity evolution. 0[L, 0] is the local luminosity function. We consider 
models where n{z) is of the form n(0) (1 -|- 2;)"° up to a turnover redshift zj^^.^^ followed by 
n(2;£„J[(l + z)/{l + -2^rn)]^'° up to z = 4.5. The luminosity evolution component g{z) = 



It should be emphasized that there is considerable degeneracy in the density and lu- 
minosity evolution of galaxies. While density evolution slides the luminosity function along 
the vertical axis, the latter slides it on the horizontal axis. However, current observations at 
mid- and far-infrared wavelengths detect galaxies only at the luminous end of the luminosity 
function, as a result of which the two are indistinguishable. This is shown in Figure 6. The 
figure also shows that evolving just the luminous end of the LLF i.e. L > 5 x 10^*^ Lq 
similar to the model proposed by Dole et al. (2000) would result in the same degeneracy 
with observations although it would result in a somewhat unphysical break in the luminosity 
function. 

There is an additional degeneracy induced by the fraction of the luminosity function 
that is evolving. Redshift measurements of ISOCAM 15 /xm sources in the HDFN-I-FF 
indicate that the majority of them are LIGs and ULIGs (Elbaz et al. 2001). Interestingly, 
when the local 60 fim luminosity function is compared to the Schechter function commonly 
used to represent the local luminosity function at visible wavelengths, an excess of galaxies 
is seen in the 60 fim LLF beyond Lgo fim > 10^^ Lq since the Schechter function drops 
faster at the bright end. On one hand it seems likely that just this excess of galaxies, 
most of which show morphological signatures of merger activity, could be evolving at high 
redshift. Alternatively, it is possible that a luminosity dependent fraction which approaches 
100% at Leo nm > 10^^ Lq of the LLF could be evolving. Unfortunately, the observational 
constraints on the faint end of the IRLF are limited since these galaxies are undetected at 
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Fig. 5. — Various 12 fim and 60 fim local luminosity functions for Ho=75 km s~^ Mpc~^. 
In the upper plot, black lines and symbols represent the 12 /im LLF from Rush, Malkan & 
Spinoglio (1993). The solid black line is the total LLF, black dash-dot line is the non-seyfert 
contribution, dotted line is from Seyfert Is, dashed line is from Seyfert 2s. The green line 
and symbols are the 12 fim LLF from Fang et al. (1998). The purple line and triangles 
are the 15 /im LLF from Xu et al. (1998) and Xu (2000). The red line is this 15 /im LLF 
converted to 12 fim based on a fc-correction derived from ISOCAM observations of 44 nearby 
galaxies. The lower plot shows the 60 /im LLF from Soifer et al. (1987) as the broken red 
line, the 60 /zm LLF of Saunders et al. (1990) as the black triangles, the 12 fim LLF of Fang 
et al. (1998) converted to 60 /xm using a linear 12 fim to 60 /xm correlation based on IRAS 
BGS data (green line) and the Xu et al. (1998) 15 /im LLF converted to 60 /im(purple line). 
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mid- and far-infrared wavelengths at z > 0.5. The correlation between mid-infrared and 
visible wavelengths being poor, the counts of galaxies at visible/near-infrared wavelengths 
cannot be used to constrain the distribution. However, we will investigate in a future paper 
if the relationship between the FIR/UV flux ratio and UV slope can constrain the evolution 
of the faint end of the IRLF. 

The principal observational constraints on the evolution of the bright end of the lumi- 
nosity function then are: 

• The ISOCAM differential number counts at 15 /xm especially the "knee" in the counts 
slope at 0.4 mJy (Elbaz et al. 1999), 

• The 15 //m extragalactic background light which has a lower limit from ISO counts 
and an upper limit from gamma-ray observations of a TeV flare in Markarian 501, 

• The redshift distribution of 15 iim sources in the HDFN-I-FF which is somewhat peaked 
at z ~ 0.8 (Aussel et al. 2001) and indicates that 45% of galaxies with 0.1 < S'15 < 
0.4 mJy arc LIGs and 20% are ULIGs with the remaining being normal and low 
luminosity starburst galaxies (Elbaz et al. 2001), 

• The spectrum of the cosmic far- infrared background between 100 and 850 iim as mea- 
sured by DIRBE and FIRAS on COBE, 

• ISOPHOT 90 yum and 170 /xm counts and SCUBA 850 fj,m number counts (Efstathiou 
et al. 2000; Hughes et al. 1998; Blain et al. 1999a; Barger, Cowie, & Sanders 1999; 
Eales et al. 2000; Dole et al. 2001). 

The /c— correction of galaxies illustrated in Figure 7 clearly illustrates the range of red- 
shifts that can be studied by observations at these wavelengths. The 850 //m observations 
are typically confusion limited at 2 mJy. However, deeper lensed surveys or using high reso- 
lution radio interferometric data can push the detection threshold down to 0.5 mJy, past the 
confusion limit. This allows detection of objects with Lj^ > 5 x 10^^ Lq out to 2; ~ 5 which 
transforms to a star- formation rate of >85 MQ/yr. On the other hand, the ISO mid- and 
far-infrared observations are typically dominated by galaxies at 2; < 1 and so can constrain 
the low redshift turnover in the luminosity function evolution. 

Figure 8 illustrates the nature of the counts if the 15 /^m luminosity function remained 
equal to the local one at all redshifts, i.e. no evolution. The first plot shows the ISOCAM 
15 iim differential counts from Elbaz et al. (1999), which include the IRAS 12 iim counts 
converted to 15 iim by Xu (2000) and the ELAIS 15 iim counts of Serjeant et al. (2000) 
renormahzed as in Genzel & Cesarsky (2000). The remaining plots show the ISOPHOT 
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Fig. 6. — The figure shows the degeneracy between luminosity and density evolution. The 
triangles are the extrapolated 60 /xm LLF of Soifer et al. (1987), the solid line is the resultant 
luminosity function at z ~ 1 assuming a density evolution of (l + z)^ while the broken line is 
the luminosity function at 2; ~ 1 assuming a luminosity evolution of (1 + z)'^. Also shown is 
the typical luminosity range of galaxies that have been detected by current long-wavelength 
surveys. 
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Fig. 7. — Predictions of flux densities at different wavelengths as a function of redshift for 
objects with an infrared himinosity of 10^^ (dash-dot-dot-dot), IQ-*^^ (dotted), 10^^ (dashed) 
and 3x10^'' Lq (dash-dot). The sohd horizontal line is the sensitivity of the deepest unlcnscd 
observation performed. The observations of lensing clusters at 15 //m are sensitive down to 
50 II Jy while those at 850 /^m are sensitive down to 0.5 mJy. 
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FIRBACK 170 //m counts, ISOPHOT, IRAS BGS and PSCz 90 //m counts and SCUBA 
850 yum integral counts. Also shown are the relative contributions from ULIGs, LIGs and 
Ljji < 10^^ Lq galaxies to the counts at different wavelengths^. Clearly, some form of redshift 
evolution in the luminosity function is required to fit the counts. It should be emphasized 
that evolutionary models should be fit to the differential counts at different wavelengths since 
integral counts tend to smooth over any subtle changes in the galaxy counts slope. This 
requires that the calibration of the data from different surveys using the same instrument 
be consistent and accurate. 

The ISOPHOT 90 fim data arc known to suffer from large calibration uncertainties 
which are not reflected in the error bars. We show in Section 3.3 that our models fit the 
170 /im counts but consistently overpredict the 90 /im counts although both wavelengths 
probe similar populations of galaxies at z < 1.2. We find that an upward correction of 
the 90 /im fiux densities by 30%, which is well within the calibration uncertainties, leads 
to excellent agreement between our models and the data. In addition, the faint end PSCz 
counts are known to suffer from incompleteness (Efstathiou et al. 2000). 

Any evolution of the luminosity function must have a turnover at some redshift Zfum 
to avoid overproducing the CIRB. In theory, the turnover redshifts z^^.^ {zturn)^ aoicxL) 
value for the slope of the density (luminosity) evolution at 2; < Ztum and the Pd{Pl) value 
for z > Ztum can be different implying that there are effectively six parameters. We however 
consider models with -2^„j= -z^rn since it is unclear why the turnover for luminosity and 
density evolution, both of which are probably induced by galaxy interactions, should be 
different. The range of values for a, {3 and Ztum selected for our models were 1 < a < 6, 
— 3 < /3 < and 0.6 < Ztum < 1-5 respectively. Evolutionary models with pure luminosity 
evolution and pure density evolution are also considered. 



3.1. Constraints on Pure Density Evolution 

We find that pure density evolution of the entire luminosity function cannot reproduce 
the counts at all the wavelengths. In this scenario, the 15 /xm counts are dominated by 
the normal/SB galaxies, not by LIGs and ULIGs, which is inconsistent with observations in 
the HDFN+FF. Secondly, the normal/SB galaxies are unable to reproduce the break in the 
15 //m counts seen at 0.4 mJy but instead produce a sharp break only at S,j < 0.2 mJy. 

However, density evolution models which evolve just a fraction of the 15 /xm luminosity 



^From now on we refer to the Lju < 10^^ Lq galaxies as normal/SB galaxies 



- 21 - 





Fig. 8. — Results for a no-evolution model with thin solid line representing the contribution 
from all galaxies, dotted line indicating ULIGs only, dashed line is LIGs only while dot- 
dashed line is for normal/starburst galaxies i.e. Lj^ < 10^^ Lq (This convention will be used 
for the remaining figures as well). Clockwise from top left 1. ISOCAM differential counts at 
15 iim. Also shown are the bright end IRAS counts of Xu (2000) which potentially have a 
factor of 2 uncertainty associated with them, 2. SCUBA 850 /xm integral counts, 3. IRAS 
and ISOPHOT 90 /xm differential counts and 4. ISOPHOT 170 //m differential counts (See 
text for references). 
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function with the fraction being <5% at Li^^m ~ 10^ Lq and approaching 100% at 15 /im 
luminosities greater than 8x10^° Lq provide reasonable fits to the data (red broken line 
in Figure 14). The best fit density evolution parameters then are an = 12.0 ± 0.5 up to 
^turn ~ ^-"^ ±0.1 followed by —0.5 < /9d < (Figure 9). The first plot shows the spectrum 
of the cosmic infrared background with lower limits from integrated counts of galaxies in 
the optical/UV from Madau & Pozzetti (2000), measurements in the near- and far- infrared 
using the DIRBE instrument (Hauser et al. 1998; Finkbeiner et al. 2000; Gorjian, Wright 
& Chary 2000; Wright 2001), an estimate of the far-infrared background from FIRAS 
(Lagache et al. 1999), and lower limits in the mid-infrared, far- infrared and submillimeter 
from counts of individual galaxies (Elbaz et al. 1999; Blain et al. 1999a; Matsuhara et al. 
2000). Also shown is the upper limit on the CIRB from TcV observations of Mrk501 (Stanev 
& Franceschini 1998) and the cosmic microwave background at A > 300 /xm. The remaining 
plots show the counts for this evolution model. 

Models with an > 13.0 rcsTilt in a significant overprcdiction of the 170 //m counts while 
«£) < 11.0 underpredicts the submillimeter counts. Changing the turnover redshift ^^^^n 
to high redshift shifts the knee in the 15 /xm differential counts to fainter flux levels and 
vice- versa. The slope of the evolution at zj^^^ > 0.7 is mainly constrained by the spectrum 
of the CIRB at A > 200 fim and the 850 fira counts. 

The 170 fim, 90 fim and 15 /xm observations all probe the population of galaxies at 
z < 1.2. However, in this density evolution model, the 15 fim counts arc dominated (~90%) 
by LIGs and have ~5% contribution from ULIGs. In comparison, the 170 fim counts are 
dominated by ULIGs at redshifts between 0.5—1 while the 90 /xm counts have roughly equal 
contributions from all three populations of galaxies. The 850 /xm galaxies at flux densities 
larger than 1 mJy are mainly ULIGs at z > 1 while LIGs between redshifts of 0.6 — 2.0 
dominate the counts at fainter flux levels. 

The pure density evolution model shown appears to overpredict the contribution from 
LIGs (~90%) to the 15 iim number counts and underestimates the bright end 850 iim counts. 
Secondly, many LIGs and ULIGs have been morphologically associated with disturbed sys- 
tems. So a steep evolution in the density of objects should reflect in an increase of the 
merger fraction which is defined as the fraction of galaxies in close pairs. Observationally, 
the merger fraction when averaged over all galaxies appears to evolve much slower with red- 
shift, approximately as (1 + zY (e.g. Le Fevre et al. 2000). It is possible though that the 
LIGs and ULIGs have a merger fraction that increases much more rapidly than {1 + zY but 
this has not been estimated since there are no clear observational signatures of LIGs and 
ULIGs at visible wavelengths. Although there is no strong observational evidence in favor of 
pure density evolution of a fraction of the luminosity function, we show in Figure 14 that it 
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does predict the same number density of infrared luminous galaxies at high redshift as other 
models and so cannot be entirely ruled out. 

3.2. Constraints on Pure Luminosity Evolution 

We have shown above that some form of luminosity evolution is required to avoid deriv- 
ing large values for the slope of the density evolution. The slope of the luminosity evolution 
a Lis strongly constrained by the mid-infrared number counts. Values of q;l < 4.5 are unable 
to reproduce the break of the differential counts seen at a 15 fim flux density of 0.4 mJy 
while ai > 5.5 results in an overproduction of the counts at bright flux densities (Figure 10). 
If the entire LLF is evolved (solid black line in Figure 14), = 5.0 up to 2;^„j=0.8 followed 
by —0.5 < Pl < 0.0, the 90 /xm differential counts are overproduced but all the other counts 
are reproduced very well. 

This is only partially consistent with the results of Xu (2000) who suggested that the 
mid-infrared counts can be modelled by evolving the entire LLF hy L{z) cx (1 + zY'^ for 
z < 1.5 and by L{z) = L{0) x 2.5^'^ for higher redshifts. We find that a luminosity evolution 
of (1 -l-^;)'^'^ up to 2; ~ 1 followed by L{z) = L(0) x 2.0^'^ overpredicts the CIRB and provides 
only marginal fits to the mid-infrared counts at the faint end. Extending this evolution up 
to z^j.^ ~ 1.5 severely overpredicts the CIRB as well as the observed faint end 15 iim and 
850 /im counts. 

The main problem with pure luminosity evolution is that the 90 jim counts at 5*1, < 6 Jy 
are overproduced but as mentioned earlier, this can be resolved by rescahng the ISOPHOT 
90 /im fiux densities upwards by 30% which is within the calibration uncertainties of the 
instrument. In addition, the break in the 15 /im differential counts in this evolution model 
is not as sharp as observed. 

As in the pure density evolution model, the counts at 15 /im, 90 /im and 170 /im all trace 
galaxies at 2; < 1.2 but the relative contributions from LIGs, ULIGs and normal/SB galaxies 
differ, with the 90 /im counts having roughly equal contributions from all three populations, 

the 170 fim counts being dominated by ULIGs and the 15 /xm counts being dominated by 
LIGs and low redshift normal/SB galaxies. Pure luminosity evolution predicts that the 
contribution from normal/SB galaxies to the 15 /xm counts between 0.1 and 0.4 mJy is 35%, 
similar to that observed in the HDFN+FF. However, the contribution from ULIGs in the 
same fiux density range is found to be only 6% in this model; a factor of 3 — 4 smaller than 
that observed. 
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Fig. 9. — Results for the best fitting pure density evolution model with only the most lumi- 
nous end of the LLF evolving (Red dashed line in Figure 14). The evolutionary parameters 
are ctu = 12.0, /5d = 0, zj^,.^ = 0.7. The upper left plot shows the spectrum of the cosmic 
infrared background: triangles represent the DIRBE/COBE results, squares with upward 
arrows are lower limits from integrated counts at different wavelengths, the dark band is the 
FIRAS constraint, hatched region represents upper limits from TeV gamma ray observations 
of Mrk501, dashed line represents the cosmic microwave background, thick solid line is the 
prediction from the model. The remaining plots are similar to Figure 8. 
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Fig. 10. — Results for the best fitting pure luminosity evolution model of the entire LLF 
(Solid black line in Figure 14) with q;l = 5.0, I3l = 0, = 0.8. (3l = —0.25 reduces 

the CIRB at A > 200 /im but slightly lowers the bright end SCUBA counts at flux densities 
greater than 2 mJy. 
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3.3. Combination of Luminosity and Density Evolution 

We have already illustrated the degeneracy and problems with pure luminosity and 
pure density evolution in the fits to the number counts. As illustrated above, both of them 
provide reasonable fits to the spectrum of the CIRB and to the differential counts at three 
wavelengths. An additional degeneracy is introduced when using a combination of luminosity 
and density evolution with only a fraction of the LLF evolving. 

The fraction of the LLF that is evolving is defined as the "dusty starburst" population. 
We consider these to be galaxies with Lb/ Li^ < 0.5. This corresponds to Lju > 10^°-^ Lq. 
All galaxies with Ljr > 10^°'^ L© are then evolved while only about 5% of the galaxies 
with LiR < 10^°"^ Lq ("normal" galaxy population) is evolved, with a smooth transition 
between the two (purple broken line in Figure 14). We also considered models where the 
dusty starburst population is defined at a larger minimum luminosity (I//r ~ 10^^ Lq) but 
were unable to find evolution parameters that could reasonably reproduce all the data. 

Our best fit model using a combination of both density and luminosity evolution for 
a fraction of the LLF as defined above is shown in Figure 11. Almost all the counts are 
reproduced quite well, with the exception of the 90 fim counts from ISOPHOT and PSCz 
which all our models consistently over estimate. We interpret this to be due to a calibration 
error in the ISOPHOT data as a result of which the published flux density values are low by 
~30%. Also shown in the figure is the surface density of galaxies per redshift bin contributing 
to the counts at different wavelengths, as derived from the model, and the observed redshift 
distribution of the ISOCAM 15 /xm galaxies in the HDFN+FF (Aussel et al. 2001). The 
size of each redshift bin is 0.2. 

In the preceding three subsections, we have shown a range of models that evolve the local 
15 /im luminosity function and fit the observed counts at mid- and far-infrared wavelengths 
as well as the spectrum of the CIRB. Ultradeep SIRTF observations at 24 //m can potentially 
break the degeneracy in these models if the counts can be determined to an accuracy of 20% 
or better (Figure 12). The range of integral source counts that we predict based on our three 
evolutionary models (density, luminosity and combination of both) are 4.7, 3.8 and 3.7 per 
for 5*1^ > 120 yuJy and 19.8, 19.3 and 11.8 per □' for S,y > 22 /iJy respectively. However, 
the integral counts can be as low as 9.1/n' at a flux density limit of 22 /iJy for models which 
are at the extreme lower limit of the uncertainty in current observations. 
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Fig. 11. — Results for a model with both luminosity and density evolution that evolves 5% of 
the luminosity function at L/^j < 10^'^'^ Lq and 100% at higher luminosities (purple broken 
line in Figure 14). ctL = 4.5, Q;i:»=1.5, /3l = 0, = —0.4, 2;4r„ = -z^^n = 
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Fig. 12. — Prediction for integral source counts seen by SIRTF at 24 /xm for the three different 
evolutionary models described in the text. Red line is for pure density evolution, black line 
is for pure luminosity evolution, solid purple line is for density+luminosity evolution while 
the broken lines are contributions from ULIGs, LIGs and Liji < 10^^ Lq galaxies to the 
counts from the density+luminosity evolution using the convention defined in Figure 8. 
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4. The Origin of the CIRB 

4.1. Nature of the GalcLxies Contributing to the CIRB 

In our evolution models, we have assumed that the contribution from AGN to the counts 
and the cosmic background is insignificant. Other evolutionary models, which assumed an 
AGN component, arrived at the same conclusion (Malkan & Stecker 1998; Rowan- Robinson 
2000; Xu et al. 2000; Franceschini et al. 2001). Observational evidence for this assumption 
comes from deep Chandra observations of the HDFN proper (Brandt et al. 2001) which 
detected 8 of the ISOCAM 15 /xm sources. However, only one of these is an AGN at z ~1 
and this was already known as such from observations at visible wavelengths (see discussion 
in Elbaz et al. 2001; Aussel et al. 2001). It should be noted that it is insufficient to have 
an AGN in a galaxy to violate this assumption but that the integrated infrared light of the 
galaxy must be dominated by an AGN rather than by star formation. However, since the 
contribution of dust obscured AGN at high redshift, beyond ISOCAM detection thresholds, is 
unknown, our results are subject to this uncertainty. A large (> 20%) contribution from AGN 
to the source counts or the CIRB will imply a weaker redshift evolution of the luminosity 
function. 

The evolution parameters in our model are constrained strongly by the ISOCAM counts 

ai z < 1.2 and by the SCUBA counts at z ~ 1 — 3. In addition, the ISOCAM counts are 
dominated by LlGs while the SCUBA counts are dominated by ULlGs. Thus, barring a 
dramatic change in the ratio of LlGs to ULlGs between a redshift of 1 and 2, we conclude 
that our models have robustly determined the evolution of the luminous end {Liji > 10^^ Lq) 
of the LLP up to 2; ~ 2. At 2; >> 2, the best constraint comes from the spectrum of the 
CIRB. Since all our models which are almost flat beyond z ~ 2 provide values for the CIRB 
that are at the upper fimit of the values observed by FIRAS, we conclude that these models 
place a strong upper limit on the estimate of dust enshrouded star formation at 2; > 2. 

Our models indicate that about 80% of the 140 /xm CIRB is produced a.t z < 1.5. 
In comparison, 90% of the 15 //m EBL, 65% of the 240 //m background and only about 
30% of the 850 //m background, is produced within this redshift range (Table 1). We also 
derive that ISOCAM galaxies brighter than 0.1 mJy contribute 14.5 nW m~^ sr~^ to the 
140 /xm EBL, while galaxies brighter than 0.05 mJy produce 16.8 nW m~^ sr""*^. This 
accounts for about 75% of the total far-infrared background (Table 2). In comparison, we 
find that the contribution from SCUBA detected galaxies brighter than the confusion limit 
of 2 mJy at 850 fim is only 3.4 nW m~^ sr^^ at 140 fim while galaxies brighter than 0.5 mJy 
produce 16.4 nW m~^ sr~^. This is because at 850 /xm fiux densities fainter than 1 mJy, the 
LIGs which dominate the ISOCAM counts and produce the majority of the 140 //m EBL 
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Table 1. Origin of the CIRB 



Wavelength 






Contribution from 2; < 1.5 galaxies 


/im 


nW m"^ 


sr 






Observed 


Model 




15 /im 


2.4±0.5 


3.2 


90% 


24 fjLm 




4.2 


83% 


140 yum 


25±7 


23.1 


82% 


240 12m 


14±3 


15.1 


67% 


850 /im 


0.5±0.2 


0.63 


28% 



Table 2. Origin of the 140 jim. EBL as derived from the models 



Source Type Contribution 



ULIGs 


25% 


LIGs 


60% 


Li [I < 10^^ Lq galaxies 


15% 


ISOCAM galaxies with 5'^(15 iim) > 0.1 mJy 


63% 


ISOCAM galaxies with 5',.(15 /xm) > 0.05 mJy 


73% 


SCUBA galaxies with 5',,(850 fim) > 2 mJy 


15% 


SCUBA galaxies with 5^,(850 /xm) > 0.5 mJy 


71% 
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Fig. 13. — Plot showing the relative contribution of LlGs (dash), ULlGs (dots) and nor- 
mal/star burst galaxies (dot-dash) to the 140 /im extragalactic background hght as a function 
of redshift. Each redshift bin is 0.2. 

contribute significantly to the SCUBA counts. The total EBL at 15 /im from our model is 
3.2 nW sr^^. The EBL obtained by integrating the observed ISOCAM counts above 
50 /iJy is 2.4±0.5 nW m'^ sr"^ (Elbaz et al. 2001), hence as much as 73±15% of the 15 /zm 
background might have already been resolved by ISOCAM. The models indicate that ULIGs 
contribute 15% of the 15 /xm EBL observed by ISOCAM above 0.1 mJy, LIGs about 65% and 
normal and low luminosity starburst galaxies the balance. In comparison, at 140 /xm, we find 
that ULIGs contribute 25% of the CIRB, LIGS contribute 60% and normal/SB galaxies the 
balance (Figure 13). Thus, infrared luminous galaxies which appear to be indistinguishable 
from normal galaxies in terms of their optical/near- infrared luminosity and form a negligible 
part of the energy budget in the local universe, dominate the star formation and therefore 
the energy budget at redshifts 2; ~ 1 — 3. 
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4.2. Infrared Luminosity Function and the LBG Connection 

The first panel of Figure 14 sliows tlie 15 fim local luminosity function (LLF) along with 
the fraction of galaxies that are evolved in the different models. For the density evolution 
model, the broken red line is evolved. For the luminosity evolution model, the whole 15 fim 
LLF shown as the solid black hne, is evolved. For the model with density+luminosity evolu- 
tion, the broken purple line is evolved. In the pure density evolution and density+luminosity 
evolution model, there is a non-evolving component with a constant comoving density which 
corresponds to the difference between the total LLF and the evolving component. 

Although the fraction of the luminosity function that is evolving and the evolutionary 
parameters are significantly different in our three evolutionary scenarios, we find that our 
models predict similar comoving number densities of infrared luminous galaxies at high red- 
shift. This is illustrated in two panels of Figure 14 which show the derived 15 /im luminosity 
function at redshifts of 0.4 and 0.8. The models also provide strong evidence for a change 
in the shape of the infrared himinosity function. The comoving number density of luminous 
infrared galaxies has to increase by more than two orders of magnitude between redshifts 
of and 1, to fit the ISOCAM and SCUBA counts. The faint end of the LLF cannot be 
enhanced by the same factor since this would lead to an over-production of the CIRB al- 
though these galaxies would be below the sensitivity limit of the long wavelength surveys. 
Lastly, as much as 85% of the far-infrared background can be attributed to infrared lumi- 
nous galaxies. This implies that the contribution from normal and low-luminosity starburst 
galaxies {Ljf> < 10^^ Lq) to the dust enshrouded star formation rate is relatively small. So, 
estimates of the total SFR by applying a constant extinction correction to all optical/UV 
selected galaxies are incorrect. We conclude that long-wavelength surveys between 15 fim 
and 850 /^m that probe galaxies at the luminous end of the infrared luminosity function 
provide a very effective way of tracing the bulk of the dust obscured star formation. 

The connection between infrared luminous galaxies and the Lyman-break galaxy (LBG) 
population is intriguing. Figure 14 shows a comparison between the LBG 60 fim luminosity 
function at 2; ~ 3 of Adelberger & Steidel (2000, hereafter ASOO), which was derived based 
on an extinction correction to optical/UV data as a function of the UV slope of individual 
galaxies, and our equivalent 60 /im luminosity function at 2; ~ 3, which we have argued 
earlier is only a strong upper limit. The agreement is extremely good considering that they 
were estimated in completely independent ways. The ASOO luminosity function predicts 
almost the same luminosity function as our estimate from the pure luminosity evolution 
model to within 50%. It is discrepant with the luminosity functions from our other two 
models by as much as an order of magnitude at the faint end but only a factor of 2 at 
the bright end. As mentioned earlier, the long wavelength surveys which constrain our 
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models are mainly sensitive to the evolution of galaxies at the bright end of the luminosity 
function at z < 2. We are unable to constrain with much certainty the evolution of the faint 
end of the luminosity function although we do place an upper limit based on the observed 
intensity of the CIRB. Furthermore, our estimates at 2; ~ 3 are only a strong upper limit 
to the number density of infrared luminous galaxies since any further evolution at high 
redshift overproduces the CIRB at A > 200 /xm while a decay in the evolution at 2; > 2 as 
[1 + {z — 2)y^ is marginally consistent with both the submiUimeter counts and the CIRB 
spectrum. Thus, we conclude that optical/UV surveys that trace the LBG population at 
z > 3, after an extinction correction factor that spans the range 2—100, provide a good 
estimate of dust obscured star- formation at high redshift. They complement the results of 
future mid- and far-infrared surveys with SIRTF which will be able to directly observe the 
dust emission of LBGs with LI G- type infrared luminosities up to 2; ~ 2.5. 

There is observational evidence for the LBG population to be distinct from the bright 
{S,y{850 /xm) > 6 mJy) submiUimeter galaxies (Barger, Cowie, & Richards 2000; Chapman 
et al. 2000). This is because most of the bright 850 fim galaxies are extreme ULIGs 
with Ljf> > 10^^'^ Lq. The ASOO observations detect only 2 out of 831 galaxies above this 
detection threshold and only 27 of their LBG sample would have a submiUimeter detection 
above the level of ~1 mJy Secondly, at 2; ~ 3, 850 /im observations would probe rest frame 
~200 /im emission. For a given far-infrared luminosity, the 60 /im luminosity shows a factor 
of 2—3 less scatter than the 850 /im luminosity among local galaxies (Figure 3) suggesting 
that the A > 200 /xm spectral shape of galaxies might potentially have a larger scatter which 
would lead to uncertain flux estimates on a galaxy by galaxy basis. Lastly, we do not find 
any scenario where the contribution from ULIGs is greater than 30% of the comoving star 
formation rate at 2; ~ 3. Naturally, the contribution from extreme ULIGs traced by the 
bright submiUimeter galaxies is even smaller. Thus, although the contribution to the star 
formation rate density from extreme ULIGs is missed in observations of the LBG population, 
their contribution is significant only at the level of <10% and hence they are less important 
to an estimate of the high-redshift dust obscured star formation. 



4.3. The Revised Stair Formation History of the Universe 

Having constrained the evolution of the mid-infrared and thereby the far-infrared lumi- 
nosity functions, we can derive the evolution of the dust enshrouded star formation rate with 
redshift. Using the equations listed in Section 2, our derived comoving SFR from all galaxies 
is shown in Figure 15 while the separate contribution from LIGs and ULIGs is shown in Fig- 
ure 16. Figure 15 also shows the absolute minimum and maximum range of dust obscured 
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Fig. 14. — Top left: The local 15 /iin luminosity function (LLF; solid line). For pure 
luminosity evolution, we evolved the whole LLF, the purple dashed line is the starburst 
component which is evolved in the density + luminosity evolution model, the red dashed line 
is the luminous component that is evolved in pure density evolution models. Top right: Total 
(evolving+non-evolving) 15 /im luminosity function for the three models at z = 0.4, Bottom 
Left: Total 15 fim luminosity function at z = 0.8, Bottom Right: 15 /im luminosity function 
at z = 3 converted to 60 /im and compared with the 60 /xm LBG luminosity function (purple 
diamonds) of Adelberger & Steidel (2000). 
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star formation rate values. The maximum values are derived from models which marginally 
overproduce the CIRB and the counts. The minimum values are derived by using an evolu- 
tionary model which is marginally consistent with the observations (Section 4.2) and by only 
considering the contribution from LlGs and ULlGs since those arc the only galaxies which 
arc directly observed at high redshift in the 15 fim and 850 surveys. Also shown is the 
SFR as inferred from direct observations at visible/UV/near-infrared wavelengths (Lilly et 
al. 1996; Madau et al. 1996; Connolly et al. 1997; Cowie, Songaila & Barger 1999; Steidel 
et al. 1999; Yan et al. 1999), SFR estimates obtained from extinction corrections to these 
observations (Madau, Pozzetti, & Dickinson 1998; Meurer et al. 1999; Steidel et al. 1999; 
Thompson, Weymann Sz Storrie-Lombardi 2001), and SFR derived from ISOCAM observa- 
tions of the Canada France Redshift Survey (CFRS) field (Flores et al. 1999). In addition, 
lower limits to the SFR from radio measurements and two points representing the complete- 
ness corrected submillimeter observations are also shown (Barger, Cowie, & Richards 2000). 
The la uncertainty in our derived rate is about 50% and is primarily dependent on the 
transformation from 15 /xm to infrared luminosities which as derived earlier has a la of 40%, 
and the transformation from infrared luminosities to star formation rate which assumes a 
Salpeter initial mass function (IMF) and also has an uncertainty of about 30% (Kennicutt 
1998). 

The dust enshrouded SFR density peaks at a redshift of 0.8±0.1 with a value of 
0.25to[\^MQyr~'^Mpc~^. In a (0.3,0.7,75) cosmology this corresponds to 6.2 Gyr after the 
Big Bang. The dusty SFR then remains almost constant up to z ~ 2 which corresponds 
to an age of 3 Gyr beyond which this value provides a strong upper limit to the amount 

of dust obscuration. This is similar to the shape of the star formation history preferred by 
Sadat, Guidcrdoni & Silk (2001) in their analysis of the CIRB. Our values arc a factor of 
2 larger than estimates at 2; ~ 1 from Ha observations by Yan et al. (1999) and a factor of 
3-7 larger than extinction uncorrected optical/UV observations at z < 2. The models are 
in excellent agreement with the submillimeter data corrected for incompleteness (Barger, 
Cowie, & Richards 2000) but are higher than the Steidel et al. (1999) extinction corrected 
points. The values we derive are systematically higher than those in Gispert, Lagache & 
Puget (2000) but within the uncertainties especially if the difference in the L/^^-to-SFR 
calibration coefficient is factored in. Our models also indicate a faster evolution at 2; < 1 
than the models of Blain et al. (1999b) which is not surprising since they did not use the 
ISOCAM data to constrain their low redshift evolution. However, our high redshift plateau 
is similar to their 'Anvil-10' model. 

Recently, Xu et al. (2000) has developed a multiparameter model where the 25 /xm 
luminosity function of Shupe et al. (1998) is partitioned into three components - starburst, 
late-type galaxies and AGN and each component is evolved independently of the other. 
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Fig. 15. — The upper plot shows the absolute maximum and minimum range of values 
derived from our model for the obscured star formation rate density in comparison with 
observed optical/UV points in a Ho=50 km s~^ Mpc~^ and go=0.5 cosmology for comparison 
to other works. Data points are from Lilly et al. (1996); Madau et al. (1996); Connolly 
et al. (1997); Cowie, Songaila & Barger (1999); Steidel et al. (1999); Yan et al. (1999) 
(black plus, blue diamonds, green crosses, inverted triangles, red crosses and purple triangle 
respectively). The dotted black line is the model of Xu et al. (2000). The lower plot shows 
our three models for the obscured SFR with the observed UV points as dotted symbols and 
the extinction corrected estimates from Madau, Pozzetti, & Dickinson (1998); Meurer et 
al. (1999); Steidel et al. (1999); Thompson, Weymann & Storrie-Lombardi (2001) as blue 
diamonds, red square, black crosses, filled red squares respectively. Also shown is the rate 
derived from ISO observations of the CFRS field (Flores et al. 1999) as filled red circles, and 
estimates from the radio and submillimeter by Barger, Cowie, & Richards (2000) as filled 
black squares. Our three evolutionary models are shown as a solid line (pure luminosity), 
broken red line (pure density) and broken purple hne (luminosity+density). We assign a la 
uncertainty of 50% to our estimate of the dust enshrouded star formation rate. We emphasize 
that our models place only a strong upper limit on the SFR at z > 2 and a drop off with 
redshift to agree with the extinction corrected optical/UV measurements is consistent with 
both the submillimeter counts and the CIRB spectrum. 
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Fig. 16. — Plots showing the contribution to the dust enshrouded star formation rate density 
from only LIGs and ULIGs as derived from our three models. A comparison with Figure 15 
shows that the contribution from normal/SB galaxies with < 10^^ Lq is relatively small 
at 2; > 0.5, ~ 5 — 30% depending on the model. Symbols are the same as those in Figure 15. 
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Specifically, they evolved the star burst population in luminosity as (1 + z)^'^ and in density 
as (1 + out to ^ = 1.5. The late-type galaxy population was evolved in luminosity as 
(1 + zY'^ while the galaxies with AGNs evolve in luminosity as (1 + z)^'^. Beyond z = 1.5, 
all the components drop-off as (1 + z)~^. Using a 25 fim to IR luminosity conversion based 
on IRAS data, wc have converted their evolution for starburst and late-type galaxies into 
a star formation rate and compared it with ours. This is shown as the black dotted line 
in the upper plot of Figure 15. We find that their derived rates between redshifts of 1 and 
2.5 are inconsistent with our models. The motivation for this peak is not clear since only 
the CIRB and the SCUBA counts place constraints on the evolution at this redshift range 
and both can be reproduced very well by an almost fiat evolutionary history at z > 0.8 (See 
Section 3). However, their evolution at z < 1 agrees reasonably well with ours since both 
are principally constrained by the 15 fim ISOCAM counts. Furthermore, their decline in the 
SFR at high redshift {z > 2) is similar but below our lower limit. 

By integrating our comoving star formation rate density over redshift and thereby cosmic 
time, we can derive the density of stars and stellar remnants and compare it with the total 
baryon density in the local Universe. Stellar lifetimes were chosen for solar metallicity stars 
(Bressan et al. 1993) while the mass of remnants was chosen using the recipe of Prantzos & 
Silk (1999, and references therein). If a Salpetcr IMF is assumed, then the model predicts 
a local density of baryons of about 1.0x10^ M© Mpc'^ which is a factor of 2 in excess of 
the value of 5±3xl0^ Mq Mpc-^ estimated by Pukugita, Hogan & Peebles (1998). The 
model also predicts that 100% of the local stars and remnants would have been produced at 
a redshift z < 2.0. If we instead use the shape of the IMP below 1 Mq suggested by Gould, 
Bahcall & Flynn (1996) which reduces the number density of low mass stars, then the 
density of stars and remnants resulting from the model is 7.5x10^ M© Mpc^^, in agreement 
with the local density. Madau & Pozzctti (2000) argued for a similar IMF based on their 
analysis of the total cxtragalactic background light. Interestingly, our model predicts that 
the local baryon density in stars and remnants, derived by integrating the star formation 
in ULIGs with a redshift distribution as shown in Figure 16, is similar to that seen in local 
spheroids, suggesting that high redshift infrared luminous galaxies may be the progenitors 
of present day spheroids. 



5. Conclusions 

A variety of observational data at mid-infrared through submillimeter wavelengths trace 
the fraction of emission from stars that is thermally reprocessed by dust. By using the counts 
of galaxies at these wavelengths, it is possible to estimate the amount of star formation that is 
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enshrouded by optically thick HII regions and thereby invisible to observations at ultraviolet 
and visible wavelengths. In addition, the spectrum of the cosmic infrared background at 
mid- and far-infrared wavelengths places an upper limit on the fraction of starlight that has 
undergone thermal reprocessing by dust. 

We have developed a set of template spectral energy distribution (SEDs) for galaxies as 
a function of infrared luminosity, which reproduce existing data at 0.44, 7. 12, 15, 25, 60, 100, 
and 850 fim from ISO, IRAS and SCUBA on nearby galaxies. The 15 fim local luminosity 
function was then evolved with redshift, considering both luminosity and density evolution 
models into account, and using the template SEDs to fit the observed counts at 15 fim, 90 fim, 
170 iJ,m and 850 /xm. A number of evolutionary models provide reasonable fits to the data 
and the spectrum of the CIRB. The principal reason for this is that all the long wavelength 
surveys are typically sensitive to only the most luminous galaxies (L/r > lO-*^^ Lq) ai z > 0.5. 
So, evolutionary models that result in similar luminosity functions at Lm > 10^^ Lq are 
degenerate. However, our models accurately constrain the comoving number density of these 
luminous galaxies as a function of redshift. In the local universe, it is these galaxies, many of 
which show morphological signatures of interaction, that show an infrared-determined star 
formation rate (SFR) that is about an order of magnitude higher than the corresponding 
UV-determined SFR. By integrating the infrared luminosity of these luminous galaxies we 
then obtain an estimate of the dust enshrouded star formation rate. The dust enshrouded 
star formation rate density appears to peak at a much lower redshift than previously thought, 
at z — 0.8 ±0.1 with a value of O.25l;o;i^-^0 yr'^Mpc'^ and remains approximately constant 
at least till 2; ~ 2. Any drop off at a lower redshift would result in an underestimate of the 
850 /xm galaxy counts. Although our models do not constrain the evolution of the faint-end 
{Lib. < 10^*^ Lq) of the luminosity function, their net contribution to the high-redshift dust 
enshrouded star-formation is negligible as can be seen in the range of evolutionary models 
considered. The evolution at 2; > 2 is constrained much more weakly. Having a constant SFR 
between redshifts of 0.8 and 4 is consistent with the CIRB spectrum and the submillimeter 
counts as is a decay by a factor of 7 between redshift 2 and ~5. However, we find that there 
is excellent agreement between our luminosity function and the infrared luminosity function 
derived from extinction correction to optical/UV observations of Lyman-break galaxies at 
2; ~ 3. This suggests that dust obscuration is significant even at 2; > 3 and that the dust 
enshrouded star formation rate is constant to within a factor of 2 between redshifts 2 and 4. 

The models also provide a census of the luminosity of galaxies that contribute to the 
counts at different wavelengths, their redshift distribution and the relative contribution to 

the cosmic infrared background at A > 5 /xm as a function of redshift. Furthermore, we 
find that ultradeep observations with SIRTF at 24 /xm down to a sensitivity of 25 /xJy 
can potentially break the degeneracy in the evolutionary models by detecting galaxies with 
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LiR ~ 10^^-^ Lq out to 2; ~ 2.5, well beyond the turnover redshift of 0.8 that is derived from 
our models. 
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